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Abstract 

We present a theory where the statistical mechanics for dilute ideal gases 

can be derived from random matrix approach. We show the connection of 

this approach with Srednicki approach which connects Berry conjecture with 

statistical mechanics. We further establish a link between Berry conjecture 

and random matrix theory, thus providing a unified edifice for quantum chaos, 

random matrix theory and statistical mechanics. In the course of arguing for 

these connections, we also observe sum rules associated with the outstanding 

counting problem in the theory of Braid groups. 
PACS numbers:05.30.-d 05.45. +b 
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1. Introduction 



To understand the theme of the paper, we present an overview of various different links 
that have been discovered in last few decades between classically chaotic systems and their 
quantal counterparts. Any study motivated to bring about this connection is what we under- 
stand here by "quantum chaos" An overwhelming number of numerical experiments 
on spectral statistics and their corresponding semiclassical analysis suggest that 
the universal features observed in chaotic quantum systems can be modelled in terms of 
random matrix theory. Apart from energy spectra, it has been found that the conjecture 
where an eigenstate of a chaotic quantum system is represented as a Gaussian random 
superposition of plane waves entails results which are found in agreement with numerical 
studies We beheve that an important step has been in establishing the result that 
this conjecture leads to momentum distribution of ideal gases, thus bringing out statistical 
mechanics 0. However, in order to bring out the puzzling results in two-dimensional statis- 
tical mechanics, it is necessary that the choice of the correlations between amplitudes of the 
eigenstates be specified. Thus, in this pursuit, we are led to random matrix theory where 
one can systematically choose the ensemble. Recently, it has been shown how one can go 
from random matrix theory to statistical mechanics - a work that has brought together 



two important statistical theories which have been, hitherto, considered quite apart. 

However, there remain many questions and we describe what we believe is an impor- 
tant one. Although one can argue for thermalization ||^ and associate a temperature by a 
suitably-defined coarse graining, is there a way we can make this association more precise 
? By the end of this paper, we hope to convince the reader that, in a certain sense, the 
way of explaining thermalization is consistent with the traditional approaches, and more 
importantly, with the second law of thermodynamics. We know that the concept of temper- 
ature associated with heat lies in a curious and subtle combination of entropy and energy. If 
we could argue that the premise that leads us to the derivation of quantum thermalization 
from the Berry conjecture also allows us to show the acceptable behaviour of entropy, we 



would have a more coherent logical scheme tied with second law of thermodynamics. As will 
be seen, this also shows at which level of description we are with respect to the projection 



methods so widely used |11 



Throughout the paper, we will be concentrating on two dimensions as that is the most 
difficult case in statistical mechanics [^^-[T^. In section 2, we give a brief discussion of 
the choice of random matrix ensemble when time-reversal and parity are broken. This is 
fundamental in dealing successfully with the problem of momentum distribution function 
and virial coefficients in section 3. The fact that quantum mechanics can be done on real 
field if the antiunitary symmetries are well-specified [llH], and, the classification theorem of 



associative division algebra |]T8| leads to three basic ensembles in random matrix theory ||17|| . 
Incorporating the violation of parity is an important step. In section 4, we unify the different 
streams of thought from quantum chaos, random matrix theory, and statistical mechanics 
by discussing entropy which is fundamental to all the three. We would like to mention that 



a recent work |jT9[ is an interesting companion of this paper. We conclude the paper with a 
summary. 

2. Random Matrix Ensemble in Two Dimensions 

In usual discussion of random matrix theory, the space dimensionality of the physical 
system plays no explicit role. Of course, it is misleading to be in that thought-frame. 
Due to the complications arising from the fact that we are working in two dimensions, we 
present here a comparative discussion about the fundamental symmetries in two and three 
(or greater) space dimensions which decisively restrict the possibilities of the random matrix 
ensemble. 

Denoting the time reversal operator by T, the position operator, q and the momentum 
operator satisfy 

TqT'^ = q, 

TpT-^ = -p. (1) 
In order to preserve the commutator between q and p, we see through 



TiT-^ = -i (2) 

(i is the square root of -1) that T is antihnear. Moreover, since 

TT^ = 1, (3) 

we say that T is antiunitary. T can always be written as a product of a unitary operator, 
U and a conjugation operator, K. On a state on apphcation of T^, we get a constant A 
times Note that 

= UKUK = UU* (4) 

which gives 

UU* = A, 

UU^ = 1. (5) 

It now follows that U = XU, and hence |Ap = 1. But then, since T is antilinear. A* = A 
which entails 

= ±1. (6) 

This then gives us, after proper introduction of angular momentum operator, the two possible 
- symmetric and antisymmetric states of a physical system which are consistent with even 
and half-odd integral spin respectively which, in turn, leads to Bose- Einstein and Fermi-Dirac 
distributions. On very general grounds thus, if a system respects time reversal symmetry, 
the Hamiltonian can be represented in terms of real or quaternion real elements depending 
on spin and rotational symmetry. If, however, time reversal is broken, the elements are 
complex, and the canonical group that preserves the Hamiltonian is unitary. Therefore, in 
three (or greater) space dimesions, a random matrix ensemble can be chosen appropriately 
satisfying the invariance under an orthogonal, a unitary or a symplectic group, and no more 

m 
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It is important to note now that the fact that we have only symmetric or antisymmetric 
states here shows that the space dimensions must be three or greater since in two dimensions, 
there is an extra phase factor under an exchange of two coordinates which leads to a fractional 
angular momentum leading to fractional statistics. In the gauge where the two particles (in 
case of a discussion of two particles one can just consider the centre of mass as one particle) 
are free, the boundary conditions get twisted. For a free charged particle of charge Q in 
a magnetic flux, the boundary condition is "^{f) ~ exp[i(integer — Q^/2Tr)(p]; we 
see that the angular momentum becomes fractional This leads to a distinction 



between clockwise and anti-clockwise rotations, which leads us to the notion of chirality and 
the associated breakdown of parity. In general, in two space dimensions, parity and time 
reversal symmetries are broken. Any choice of a random matrix ensemble must be consistent 
with this. 

Since time reversal is broken, it follows from the foregoing discussion that the Hamilto- 
nian matrix of the system will be complex, invariant under a unitary group. Breakdown of 
parity is new, however the answer is in the boundary condition. Thus, we are led to a chiral 
unitary ensemble. 

In our present context of many-body system whose all the eigenstates we cannot know 
exactly due to practical limitations (even if it is possible in some cases, we deal with the 
situation where a statistical study is the viable option), following Srednicki, we write a 
random pure state as a superposition of some basis states with amplitudes which are random. 
By the randomness of the amplitudes, we mean that they satisfy some correlation functions 
which we will write in the next section. The randomness in the amplitudes makes the pure 
states of the system also random. We have assumed that the system is isolated. 

The randomness in pure state can also be interpreted by weighting the eigenvectors 



by a measure invariant under unitary transformations, Un. By considering the unit complex 
N-sphere as a homogeneous space of Un, then again as Un itself but organized into cosets, 
the measure is seen to be the Haar measure on U at, thus unique. It is from this interpretation 
that we will discuss the entropy of the subsystem where we will note a connection between 



randomness in pure state and random matrix theory, however in that case it will be applied 
to the density operator of the subsystem which resides in the isolated system with a Hilbert 
space of lesser dimensionality. 
3. Momentum distribution 

Let us consider a system of hard spheres ('discs' in two dimensions), each of radius a, 
enclosed in a box of edge-length L + 2a. Centres of two hard spheres Xi and Xj are such that 
\xi — Xj\ > 2a. The canonical pair of coordinates describing these particles are {X,P) where 
X = {xi,X2, - ■ ■ ,Xn), P = (pi,P2, ■ ■ ■ ) PW)- Energy eigenfunctions, 'ipai^X) corresponding to 
eigenvalue Ea vanish on the boundary of the enclosure. A typical eigenfunction is irregular, 
with a Gaussian amplitude distribution and the spatial correlation function of the same is 
consistent with the conjecture of Berry which allows us to represent this eigenfunction as a 
superposition, following Srednicki 0: 

= N^j d''''PA^{P)6{P^ - 2mE^) e^^'^ (7) 

with Na given by the normalization constant, and A'^s satisfying the two-point correlation 
function 

5'^^iP-P') 

d denotes the number of coordinate-space dimensions. The average in @ is a matrix- 
ensemble (ME) average which originates from the fact that the hamiltonian, H of the system 

in consis- 



belongs to an ensemble of matrices satisfying associative division algebra ||17| , p!8 
tency with quantum mechanics. The eigenstate ensemble (EE) used in is nothing but a 
consequence of underlying matrix ensemble in RMT, the eigenfunctions then satisfy all the 
properties numerically observed and analytically represented in (|^, (H) The correla- 

tion functions (^ decide whether time-reversal symmetry is preserved (A* (P) = Aa{—P)) 
or broken (y4*(P) ^ Aa{—P)), accordingly the corresponding matrix ensemble belongs to 
OE or UE respectively. As noted in 0, the higher-order even-point correlation functions 
factorize and the odd-ones vanish. A very important aspect of the ansatz (|^), (H) is that the 



— * 

Wigner function corresponding to ipai^) is microcanonical, or, is proportional to 8{H — Ea) 
which, in a sense, incorporates ergodicity. We note here that, starting from an ansatz very 
similar to the one above, it is possible to obtain the quantum transport equation where 



it is important to relate a given quantum state with the admissible energy surface in phase 
space; thus the above ansatz is in conceptual agreement with the ergodic aspect of many- 
body system. Moreover, this choice fixes the Thomas-Fermi density of states naturally. It 
now becomes important to emphasize that we must restrict ourselves to dilute gas of hard- 
spheres and also assume that the size of sphere is much lesser than the thermal de Broglie 
wavelength. Thus, the ansatz establishes, in fact, a link between RMT and statistical me- 
chanics. We now incorporate the case of two dimensions which otherwise presents enormous 
difficulties. 

In two dimensions, the solutions of the Schrodinger equation, ip{xi,X2,- ■■ ,xn), under 
an exchange of two coordinates of particles satisfies 

■ ■ ■ , Xj, ■ ■ ■ , Xj, • ■ • , Xisf) 

= e'''"'i){xi, ■ ■ ■ ,Xj, - ■ ■ ,Xi, - ■ ■ , xn) (9) 

where u is arbitrary and defines statistics. For = and = 1, with (H), one gets the Bose- 
Einstein and Fermi-Dirac distributions. This non-trivial phase and the resulting boundary 
condition arises from the fact that the effective configuration space, has a fundamental 
group, TTi{Mjf) = [Q, the Braid group of objects which is an infinite, non-abelian 
group. Bn is generated by {N — 1) elementary moves ai, ■ ■ ■ ,aN-i satisfying the Artin 
relations, 

atai+iai = a^+iaiai+i (i = 1, 2, ■ ■ ■ , - 2) 

GjOi = OiOi, \i- 3\ > 2 (10) 

the inverse of Oi is 0"^"^, the identity is denoted by I, and the centre of B^ is generated by 
{o\02 ■ ■ ■ aN-i)^ ■ The multivaluedness of the eigenfunction originates from the phase change 
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in effecting an interchange between two coordinates xf^ and xf'^ (superscripts refering to 
components) which can be expressed as 

i<j 

(f)ij = tan-1 (11) 

The description adopted by us here is referred to as the Anyon Gauge. It is important 
to reahse that a set of coordinate configuration can be reached starting from some initial 
coordinates of particles in an infinite ways, each possibility manifested by an action of an 
element P E B^. 

The connection between initial and final sequences is given by (^, via the character 
xiP) of the specific element. Thus, to every P G B^, we can associate the affected partial 
amplitude ipaiP '■ x) p5|. With one-dimensional unitary representation of the braid group, 
the rudiments of quantum mechanics allow us to write 

$«(X) = X{P)MP ■ X) (12) 

where ipaiP '■ X) is the probability amplitude associated in changing a configuration X to 
(/3 : X) - a configuartion after the action of P on X. The wavefunction $q,(X) is to be 
understood as appropriately normalised. The ansatz for Vipa{P '■ X) is now 

V^Ij^{P ■.X) = N^j (f^PA^{P : P)5(p2 _ 2mE^)e^^-^ (13) 
with Aa{P : P) satisfying 

{^;(ft : P.)Mlh : A))m. = (14) 

^l-' 1 ~ -'2 ) 

{Pi,P2 G Bm), and Aa{P) satisfy the twisted boundary conditions, 

^a{pir--,Pi,---,Pj,---,PN) 
= e^^'Ao.ipu ■ ■ ■ ■ ■ ■ ■ ■ ■ ,PW) (15) 
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The question now is in specifying exactly what the matrix ensemble is in this case? The 
form of (|n|) with A^s not restricted to real, takes into account the T-breaking, and (|T5D 



makes the ensemble handed or chiral as a result of P-breaking. Thus (p!^)-(|T5|) gives the 
complete description and the ME is, in fact, the chiral- Gaussian Unitary Ensemble (ch-GUE) 
as discussed in the previous section from general considerations. It can be easily shown 
that the Wigner distribution is 

p2 



{p:{X,P)),,j, = n-'h-'''6{—-E, 



— 1 rriL'^Ea 

" iV!r(iV)E/ 27:h' ' ■ ^ ^ 

For the momentum distribution, we need to evaluate the ME- average of $*(P)$*(P") 
with $ = V^a- With the above ansatz and conditions supplementing it, this average is 



^(p) = (i>:,(p)i>!(p")) 



ME 



h^''5^.,N^,N, E E X*(/5i)x(/32)5(P'-2mP, 

n,m=0 /3i{m) P2{ri) 



X( 



P=P" 



(17) 



where 



4^(Q) = h-'"" f rf^TV^exp (^Q ■ X); (18) 

J Domain,T> 7l 



P is identified with P" after the sum is performed. 
With (p!7|), the momentum distribution is given by 



S dpi--- dpNJ^{P) 

which formally completes the deduction. However, an exact evaluation of this is very difficult 
and the difficulty is coming from counting of irreducible words formed by the a's. To make 
the precise connection, we give derive the result upto 0{h^ /L?'), an order that is enough for 
second virial coefficient. 

In deriving the momentum distribution, we have to consider all the exchanges that lead 
to contributions giving second virial coefficient. With N generators, we have characters e^'"" 
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and e~*'^^ leading to a combination, cos(A^7ri/). With one of the momenta fixed in the above 
integral, we look for elements of Bf^ such that two momenta are interchanged restoring all 
other momenta to their labellings. All these elements contribute upto 0{Ti^ /Li^). There are 
two kinds of terms with any um (which denotes the elements oi with M generators): 

(i) one where pi changes, 

(ii) one where p[ does not change. 

We first insert a notation which will be used in sequel, viz., the integral, 

(20) 

In case (ii) , we have typically 

Soipi -Pi)Sd{P2 -p2)--SD{Pj -Pi)SDiPi -Pj)--SD{PN -Pn)- (21) 

This leads to the value of the integral, 

2 V-2)(2m^a - Pi) (J n,^{N)x{aM) (22) 

where T^^m denotes the number of elements composed by M generators that contribute to 
0{h? / L'^), or just an interchange between two momenta but not pi, and, x{^m) denotes the 
corresponding character. 

In case (i), we have typically 

Sd{P2- Pi)Sd{Pi -P2)--5d{pn - Pn)- (23) 
which leads to the integral evaluating to 

/^x 2(iV-l) 

hiN-2){2mE^ - 2pl) (^-j QsAN)x{Bm) (24) 

where Qffui^) denotes the number of elements of Sjv composed by N generators contribut- 
ing to 0{}{^ jl?') that involve an interchange with p\. 

For each ou x{^m), we can find X*{^m)-, and (Q<tm>'^<?m) = (Q'?!,' ^'?m)- 

Thus, for fixed M, case (i) gives 
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2(7V-1) 

/2(iv-2)(2mE„ - 2pl)Qff^{N) (xI^m) + X*{^m)) ; 



(25) 



and case (ii) gives 



L 



h 



2(N-l) 



'2(Af-2) 



{2mE^-pl)ng^{N){x{aM)+X*{^M)). 



(26) 



Because only two momenta are interchanged, the total contribution of elements of 5jv formed 
by M generators is 



hj 



2(iV-l) 



E 

k=-M,-M+2,...,M-2,M 



l2iN-2){2mE^-pl) cos{jiku)nf\N) 



(27) 



for 



—iMnu —i{M—2)'Kv iMirv 



(28) 



and, TZgj^is just TZ^^\n). If we integrate over pi, we obtain the normalization factor. For 
this, we have terms that lead to exchange as discussed above, and also the operation of 
identity of where no momenta are changed. To begin with, we have the integration of 
the term without identity, and the result is 



^^2(7V-l) 



E 

k=-M,-M+2,...,M-2,M 



+272(iv-i) {2mEa) cos(7rA;z/) qI*^^ {N) 



(29) 



Denoting by Vffj^ {N) by the number of elements of Bjy composed of M generators con- 
tributing to 0{hP / LP) - the identity, integration over p2---PN gives 



-j hiN-i){'2mE^-pl)V3^{N)x{aM). 



(30) 



As above, we have &m and ct^, so this integral reduces to 



h 



2N 



^ 2h^M-i){2mE^ - pi) cos{7rku)V^{N). 

k=-M-M+2,...,M-2,M 



(31) 



To get the contribution of identity to normalization, we now integrate this over pi to obtain 
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//■n 2iV 

{-r E 2l2N{2mE^)cos{7rkiy)V^{N). (32) 

k=-M,-M+2,...,M-2,M 

For the second virial coefficient, if M = 2m (m=0,l,2,...), tfie contribution goes to tfie 
term involved in identity, and, if M = 2m + 1 (m=0,l,2,...), tlie contribution is 0{h?' / LF'). 
Tlie sum over elements of Sjv can be substituted by a sum over m. All put together, in the 
term which gives normalization, we have 



V E E 272iv(2m^a)cos(7rfcz/)Pf (iV); (33) 



m=0 fc=— 2m,— 2m+2,...2m 



and 



7,2 \ .rs2(Ar-l) oo 

0(- : - E E [/2(iv-i)(2m£;«)cos(7rH<W 

' ^'''^ m=0 fc=-2m-l,-2m+l,...2m+l 

+/2(iv-i) (2m£;,) cos(7rA;i/) Qf (TV)] . (34) 



For the numerator of (19), with one momentum pi fixed and integrating with respect to all 
other momenta, we get the following results : 

.rs 2Ar oo 

Mr) E E 272(iv-i)(2m£;,-p?)cos(7rHn''W; (35) 

m=0 fc=-2m,-2m+2,...2m 

and 



O 7^ Ml- E E [/2(iv-2)(2m£;„-pDcos(7rH7ef(iV) 

yiv y k=-2m-l,-2m+l,...2m+l 

+2/2(Af_2)(2mE„ - 2p2) cos(7rA;t/)Qf (iV)l . (36) 



For large N, 

/2(iv-i)(a;) ~ /2iv(a;) (37) 
as /Ar(x) is just the volume of an N-dimensional sphere of radius x. Let us define 

A 



2mkBTa 
_ 1 

^ ^ ^^^^ 
With these, the terms for the normalization factor can be re-written as : 
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0(1): -) 2hN{2mEo)Y. E cos{'Kku)V^ {N)- (39) 

^""^ m=0 k=-2m-2m+2,...2m 



and 



^ 2(iV-l) 

: (^-j l2N{2mEa)B 



E [cos(7rA;i/)7ef (7V) + cos(7rA;i/)Qf (AT)]. (40) 

m=0 fe=-2m-l,-2m+l,...2m+l 

Similarly, the terms corresponding to te numerator of (19) can be re- written as 

Oil) : - 2h^{2mEo)Be^Y>{-A)Y. E cos{7Tkiy)V^ (N); (41) 

^'^^ m=0 k=-2m,-2m+2,...2m 

and 



//,2\ .r^2(7V-l) 



E [cos(7^A;i/)7^f (A^) + 2cos(7rA;i/)Qf(Ar)]. (42) 

m=0 k=-2m-l,-2m+l,...2m+l 

Eqs. (41) and (42) combine to give the numerator of (19) which we call Ti, and, Eqs. 
(39) and (40) combine to give the denominator of (19) which we call ^2- Thus the ratio 
of to ^2 gives the momentum distribution upto 0{h? / E^). Now, after a straightforward 
arrangement of all the terms , we get 



F(pO = (27rmA;r)-^exp(-^){l 

L_(2e- - 1)G(7V, v) + 0(^)} (43) 



where 



, E^=0 Er='^-\^-i(e.en) Q^K^ W C0s(7rX^) 

GlA/, i/) = , (44) 

1 + 2 E^=i Er=-2n^(o..) cos(7ri^i.) 



Qx^^ is the number of elements in composed of 'm' generators whereby the momentum pi 
is interchanged with another momentum yielding a character expiiTrKu) (or exp(— iTrXz/) 
since Qk^\n) = Q^(A^)); pI^\n) is the number of elements in B^ contributing to 
identity with a character exp{iTrKu) ( or exp(— ivri^i/)). Temperature is introduced above 
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via the ideal gas law, Ea = NkT^. Unfortunately though, this counting problem stands open 
today |]2^. It is very important to note that the ansatz ([T3|)-(|15|) for the special case when 



af = 1 for all i where reduces to symmetric group, Sn, the well-known Fermi-Dirac and 
Bose- Einstein distributions follow. On evaluating pressure, 11 from (14), denoting area of 
the enclosure by A, we get UA/kT = 1 - {2A)-^\'^G{N,u), with = h'^{2nmkT)-\ We 
immediately see that ^(A^, 0)/(2iV) and ^(A^, 1)/(2A^) are 2^3/2 and -2^3/2 respectively 
yielding the second virial coeffiecient for the Bose and Fermi gases |2^. For the fractional 



case, with u = even number, 2j + 6 ("boson-based anyons"), comparing our result with JTS 
we get the Sum Rule : 

- 2~^/^N-'G{N, u)X^ = (-1 + A\6\ - 26^)XyA, (45) 



the right hand side belongs to [|T5[. It is important to note that our deduction is non- 
perturbative and in principle, we can get expressions for higher-order virial coefficients also 
p9| . To understand this, we observe that the relation (14) connects two momentum config- 
urations of N particles, and not just the momenta of two particles. Thus, it contains infor- 
mation that can lead to all virial coefficients. For example, for the third virial coefficient, 
we need to evaluate contributions to F{pi) when three momenta out of N are interchanged. 
The denominator of (19) contains those interchanges which braid three strands in such a 
way that the initial configuration of momenta is preserved whereas the numerator of (19) 
contains those which exchange the momentum assignment on all three strands. We have 
done the calculation and the third virial coefficient is expressible in terms of the specific 
counting problem of -Bat. Here, in order to convince the reader, it suffices to make a compar- 
ative discussion with the existing calculation. For this, we write down the total contribution 
to the momentum distribution due to a triple interchange emerging from the elements of 
Bj^i formed by M generators, 

r\2{N-2) 2M 2 

- E -hiN-3){2mE^ - f,) cos(7^A;z/)7^f (AT) 

+ 2/2(^_3)(2mE„ - 3p?) cos{Txku)St\N), (46) 
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where S^{N) {TZ^{N)) are the number of elements of that (do not) change the mo- 
mentum pi. Id{x) denotes the volume of a D-dimensional hyper-sphere of radius x. The 
reason we give this result here is to show that (46) is a Fourier series with harmonic terms 



like cos27rz/, cos47rz/, and so on, in complete agreement with the conjectured form |2^. It is 



becoming evident from the Monte Carlo calculations BO] that the third virial coefficient is 



a series with terms as sin^ vrz/, sin^ vrz/, and so on. Our formal result is thus in consonance 
with these works. Also, we mention that (45), (46) and the Monte Carlo estimates provide 
a non-trivial hint on the counting problem itself. 

4. Average entropy of a quantum subsystem - averaging over random eigen- 
states, and, over random Hamiltonians 

In this section, our discussion will not be restricted to two dimensions. Also, the subject 
will be entropy which is apparently different from the previous section. The problem that 
we address is a quantum version of the Ehrenfest urn model |^ as first considered in p2 



Consider a system AB with Hilbert space dimension mn and normalized density matrix 
p (a pure state p = \ip >< ip\ if = p) [^,^]. Now we divide this system into two 
subsystems, A and B, of dimension m and n respectively. The density matrices of A and 
B, respectively, pa and pb, are obtained by partial tracing of p over B and A respectively. 
We assume that A and B are quantally uncorrelated, i. e., p = pA ® Pb- If I'?/' > is chosen 



at random, what is the joint probability distribution of eigenvalues of pa ? Following P2 
"random" refers to unitarily invariant Haar measure which, in this case, turns out to be 
hyper- area of the unit sphere ^^mn-i^ ^^le factor 2 coming from the fact that {ip > has 
mn complex entries (or 2mn real entries). The objective is to study the average entropy 
of A, < Sa > (=-tr Pa log pa) over the probability distribution of eigenvalues (which are 



probabilities) of pa- The result of this calculation, conjectured in |^ and proved first in 
is 

mn 1 m — 1 

<Sa>= E (47) 

fc=n+l 

The random pure state can be written (for the system we are considering) as (7,8) in 

15 



three (or greater) space dimensions, or, as (13-15) in two dimensions. These equations can 
be interpreted as choosing a pure state at random for a specific choice of amphtude, A^. 

If one calculates the average of trace of p\, first, over homogeneously distributed unit 
vector in mn-dimensional Hilbert space, and then, over random Hamiltonians, the two an- 
swers are only different by one bit |^ . The values are almost the same as one corresponding 



to the answer when the entropy will be maximal, i.e., logm. This brings about the random 
matrices as all that is being done here about averaging over random Hamiltonians is what 
is done in random matrix theory. Thus the average entropy of a subsystem follows from the 
random matrix hypothesis about AB - the statement becomes exact when m <^ n. Indeed, 
the connection of random matrix theory and statistical mechanics is when the size of the 
system is large where it means then that the number of particles is large to be consistent 
with thermodynamic limit. 

It is very interesting to note that if the pure states of AB are random, the probability 
distribution of eigenvalues of pa is just the one-point correlation function corresponding 
to the (generalized) Laguerre unitary ensemble of random matrices Since the one- 

point correlation function (average level density) is the same for orthogonal, unitary, and 
symplectic ensembles the answer for the average entropy will remain the same. Let us 



remember that the average level density is, in principle, a function of the size of the matrices. 
The fact that we must discuss entropy in the context of systems in thermodynamic limit 
is what makes the entropy same for all the ensembles. Since the fluctuations on top of 
the average density become signiflcant with decreasing sizes, we expect to observe their 
interesting effect on the entropy. We now present results that prove these remarks. 

We begin by recalling that the average over all pure states of AB, in unitary Haar 



measure, of the spread of eigenvalues of pa is |22 



<a'> =<-22[Pi > 



m ^ \ m 

_2 

m 



mn + 1 

Clearly, the case n=l corresponds to the situation when pA is also pure, then we have 
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(48) 



the ratio of < a > to < cr^ax > gives the measure of "purity" of the subsystem A. Since 
we have noted above that the same answer can be obtained by averaging over the one- 
point correlation function of the Laguerre unitary ensemble of random matrices, we expect 
that apart from the leading term which is logm for the entropy as this corresponds to 
equipartition, the "defect" term must show the signature of different ensembles. To this 
end, we start by writing the entropy, 

S = -Y^Pilogpi, (50) 

i 

and Taylor expand each pi about 1/m. With 

Pi = , 51) 

m m 

we can write 



S = log m 

m 



1=1 1=1 

2 q;^^„ ^ ^2 



(52) 



which is convergent if \qi\ < 1, i.e., if < Pi < ^. Since < o" > is small for large n, most 
of the measure will lie with pi < ^. It is plausible that S = logm — defect, and that the 
defect is well approximated by 



1 ^ . 1 



< defect > = q. 



2 -^22 
— ma 



2m 2 
1 — 1 



(53) 



2 mn + 1 

We can now find the defect for the case where we integrate over orthogonally invariant 



Haar measure |^ and symplectically invariant Haar measure. The difference is that these 
correspond to and respectively. After the same steps, we get 

/ aM |m + 1 



i '^max I ^mn -\- \ 



(54) 



where (3 is the co-dimension of level crossing, respectively 1,2, and 4 for orthogonal, unitary, 
and symplectic ensembles of random matrix theory. The entropy is given by 
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(m-1) fm+1 
5 = logm -f — . (55) 

This result shows that the entropy is almost maximal for n large enough, and that the 
finite-dimensional effects show up the dependence on the global symmetries of the system. 

We wish to note that the arguments we have used are quantum mechanical. The results 
of this section show that randomness in eigenstate, which follows from random matrix theory, 
encompassing the ergodicity of the classical system leads to entropy of a subsystem which 
is maximal. To notice the differences between the chaotic quantum systems belonging to 
different random matrix ensembles, we need to study the entropy as a function of the Hilbert 
space dimension as is clear above. We conclude this section with the mention of a recent 
work where a related study is carried out on periodically kicked top ll^ . 
5. Concluding Remarks 

We have shown in this paper that an ansatz where eigenstates are written as random 
superpositions of plane waves (or some basis consistent with the boundary conditions) for 
systems whose classical analogues are chaotic is equivalent to a random matrix hypothesis. 
In an important work 0, the connection between the ansatz and statistical mechanics in 
three dimensions was brought out. Since the group that governs exchange symmetry in two 
dimensions is an infinite, nonabelian one whose special case is the permutation group (the 
exchange group in three dimensions), our treatment and results in Section 3 are generaliza- 
tions of 1^. It was possible for us to do this only because we realised what the random matrix 
ensemble should be in two dimensions. This is the reason for Section 2 where a comparative 
discussion about the relevance of space dimensions is given to the random matrix ensem- 
bles. In consistency with the expectations, the ensemble in two dimensions is chiral unitary 
ensemble. It is well-known that choosing a specific nature of randomness (e.g., Gaussian) 
gives then the average density of states which is not realistic. This can be treated with the 



Dyson Brownian motion model pO[ where any realistic density of states can be modelled. 
An interesting relation between a generalized Brownian motion model and a semiclassical 
reasoning of universality has been recently worked out pM] by a generalization of the theory 



of level dynamics. 

The fact that both the ideas, one of random pure state of an isolated system, and, that of 
this system being governed by a random Hamiltonian, give rise to the distribution functions 
and virial coefficients correctly suggest that there may be a connection between the two. In 
Section 4, we have described a way that we see most clearly (as of now) in the context of 
entropy of a quantum subsystem. We believe that the arguments developed here provide a 
common ground to seemingly different themes, leading to the second law of thermodynamics. 

As shown in Section 3, the number of words formed by M generators of the Braid group 
satisfy a sum rule which comes from our calculation of the second virial coefficient. The 
parameter, u has an analogous partner in quantum chromodynamics [|15| and we conjecture 
that the anyon gas discussed here and the z/vr-parametrised quantum chromodynamics belong 
to the same universality class of chiral Unitary Ensemble of RMT. 
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